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jj 要 : 脉冲 星 周期 跃 变 是 一 种 军 见 的 现象 ， 是 研究 其 内 部 结构 的 探 针 。 针 对 2019 年 2 
月 至 12 月 ， 国 家 授时 中 心 昊 平 观测 站 40 米 射电 望远镜 在 脉冲 星 计时 观测 中 ,监测 Crab 脉 冲 星 
的 数据 ， 采 用 脉冲 星 计 时 方法 ， 用 TEMPO2 拟 合 程序 进行 分 析 。 结 果 表 明 Crab 脉 冲 星 在 2019 
年 7 月 23 日 (MJD58687) 附近 发 生 了 一 次 周期 跃 变 现象 , 该 跃 变 自转 增 量 为 Avs= 5.33 (4) 
X10 ”Hz， 自 转变 化 量 为 Av, /v 2179 (D x10? ， 并 伴随 着 恢复 系数 O ~ 0.88 的 指数 恢 
复 过 程 。 此 次 Crab 脉 冲 星 周期 跃 变 的 监测 及 处 理 , 证 实 了 国家 授时 中 心 40 米 射电 望远镜 对 脉 
冲 星 的 监测 性 能 ， 同 时 为 研究 周期 跃 变 的 产生 机 理 积 累 了 样本 。 
关键 词 : AAKE; AFE; EREZA; THRE 
中 图 分 类 号 : P161.1 文献 标识 码 : A 文章 编号 : 


脉冲 星 是 一 种 具有 极 强 磁场 的 高 速 自转 的 中 子 星 ， 其 长 期 稳定 度 可 妇 美 原子 钟 站 。 部 分 
脉冲 星 自转 表现 出 计时 噪声 (timing noise〉 和 周期 跃 变 (glitch) 这 两 类 不 稳定 性 。 计 时 品 
声 广泛 存在 于 脉冲 星 中 ， 是 一 种 连续 的 、 小 幅度 的 波动 现象 四， 在 具有 较 大 周期 导数 的 脉冲 
星 中 表现 更 为 显著 , 大 部 分 计时 噪声 具有 准 周 期 特征 ; 周期 跃 变 则 表现 为 自转 速度 突然 增加 ， 
通常 伴随 着 时 长 几 周 到 几 年 不 等 的 恢复 过 程 趾 ， 周 期 跃 变 事件 比较 罕见 ， 是 一 种 偶发 现象 ， 
目前 无 法 进行 预知 内。 在 目前 已 知 的 2872 颗 脉 冲 星 中 ”， 共 监测 到 190 颗 脉 冲 星 的 567 次 周 
期 跃 变 ”。 

在 已 监测 到 的 大 多 数 周期 跃 变 事件 中 , 没有 发 现 辐射 特征 的 明显 变化 。 因 此 ， 周 期 跃 变 
很 可 能 与 脉冲 星 内 部 结构 的 变化 有 关 。 目 前 关于 周期 跃 变 的 产生 机 制 主 要 有 两 种 理论 模型 : 
星 震 模型 名 和 涡流 模型 四 。 星 震 模型 认为 ， 周 期 跃 变 是 脉冲 星 壳 层 形变 而 导致 的 现象 ， 即 脉 
冲 星 自转 的 稳定 减 慢 会 破坏 引力 与 离心 力 的 平衡 , 脉冲 星 发 生 突然 收缩 , 使 得 脉冲 星 自转 突 
PRIN. 涡流 模型 认为 ,周期 跃 变 是 在 脉冲 星 自转 减 慢 过 程 中 ， 内 部 超 流 体 和 外 部 壳 层 的 角 
速度 存在 差 速 旋转 ， 当 壳 层 的 旋转 速度 比 内 核 慢 且 耦 合力 不 超过 某 种 极限 值 时 ,， 镶 庶 在 外 过 
品格 中 的 涡流 只 能 向 外 运动 ， 并 把 角 动 量 传 给 外 这， 于 是 形成 脉冲 星 自转 速度 加 快 现象 。 
星 震 模型 可 能 适用 于 小 幅度 的 周期 跃 变 ， 涡 流 模型 可 能 适用 于 大 幅度 的 周期 跃 变 。 

1968 年 发 现 的 Crab 脉 冲 星 (PSR B0534+21 或 PSR J0534+2200) 形成 于 公元 1054 年 ， 是 
蟹 状 星云 超新星 产生 的 致密 星体 ， 同 时 是 目前 已 知 最 年 轻 的 脉冲 星之 一 四。Crab 脉 冲 星 的 自 
转 周 期 约 为 33ms, 周期 导数 约 为 4.2x103， 有 约 4x102 G 的 强 磁 场 , 具有 多 波段 辐射 特性 中。 
但 由 于 Crab 脉 冲 星 比较 年 轻 ， 内 部 结构 还 不 太 稳 定 ， 因 此 ,研究 Crab 胀 冲 星 的 周期 跃 变 事件 
可 作为 研究 中 子 星 内 部 结构 的 探 针 ， 具 有 重要 科学 意义 。 

本 文 将 介绍 中 国 科学 院 国 家 授时 中 心 吴平 观测 站 40 米 望远镜 监测 Crab 脉 冲 星 ， 在 2019 
年 7 月 23 日 前 后 自转 周期 发 生 跃 变 的 情况 。 第 一 章 介 绍 了 用 于 本 次 研究 的 观测 数据 ， 第 二 
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首 述 了 数据 处 理 流程 ， 第 三 章 对 计时 人 处理 结果 进行 分 析 ， 在 第 四 章 给 出 结语 。 


1 观测 


观测 使 用 中 国 科 学 院 国家 授时 中 心 National Time Service Center, NTSC) 的 40 米 射电 望 
远 镜 (Haoping Radio Telescope, HRT) 。 该 望远镜 建 于 2014 年 ， 位 于 西安 以 东 约 100 公 里 的 
洛 南 县 吴平 观测 站 ， 坐 落 在 四 面 环 山 的 秦岭 山脉 中 ， 电 磁 环 境 优良 。HRT 是 一 架 全 可 动 的 卡 
塞 格 伦 式 单口 径 望远镜 ， 由 全 面板 反射 面 组 成 50。 脉冲 星 计 时 观测 系统 配备 工 波 段 常 温 接 
收 机 和 FPGA + GPU 架构 的 数字 终端 。 其 中 ， 接 收 机 带宽 为 800MHz， 在 2019 年 11 月 之 前 为 
右 旋 圆 极 化 的 单 极 化 接收 模式 ， 此 后 升级 为 双 圆 极 化 观测 模式 。2019 年 7 月 开始 新 增 相 干 消 
色散 模式 。 系 统 噪声 温度 约 100K。 

吴平 观测 站 望远镜 对 Crab 脉冲 星 的 观测 为 2-5 天 进行 一 次 ,观测 时 长 10-75 分 钟 不 等 。 
于 本 次 Crab 脉冲 星 周期 跃 变 研究 数据 的 时 间 跨 度 为 2019 年 2 月 21 日 至 2019 年 12 月 8 
日 (MJD58534 - MJD58825), 其 中 包括 非 相 干 消 色 散 和 相干 消 色 散 中 两 种 观测 模式 的 数据 。 
在 非 相 干 消 色 散 模式 下 ， 时 间 分 辨 率 为 10.24us， 采 用 1024 个 频率 通道 和 1024 个 相位 进行 
数据 记录 ; 在 相干 消 色 散 模式 下 ， 时 间 分 辨 率 为 1.28us， 采 用 1024 个 频率 通道 和 1024 个 相 
位 进行 数据 记录 ， 子 积分 时 间 为 10s。 


2 数据 处 理 流程 


计时 观测 数据 处 理 采 用 预 处 理 软件 PSRCHIVE 和 计时 分 析 软 件 TEMPO2。.PSRCHIVE02 
是 一 款 科学 数据 分 析 的 开源 软件 , 应 用 于 脉冲 星 天 文学 研究 , 可 实现 脉冲 星 计 时 观测 数据 的 
校准 、 统 计 分 析 与 模拟 以 及 可 视 化 TEMPO203 是 为 脉冲 星 计 时 阵列 项 目 开发 的 一 款 计 时 软 
件 ， 使 用 国际 天 文 参 照 系 “ICRS) 并 遵从 IAU2000 决议 ， 用 于 拟 合 周期 跃 变 参数 ， 计 时 精 
度 可 达 纳 秒 级 别 。 
本 次 数据 处 理 流程 如 图 1 所 示 ， 对 ART 观测 所 得 折合 模式 数据 ， 首 先 使 用 PSRCHIVE 
对 数据 进行 消 干扰 ; 再 对 消 色散 后 的 子 积 分 轮廓 在 在 时 间 、 频 率 和 偏振 上 进行 个 加 ， 得 到 平 
均 脉冲 轮廓 ; 为 了 得 到 更 准确 的 脉冲 信号 到 达 时 间 ， 通常 把 所 有 观测 到 的 脉冲 轮 廊 合 加， 得 
到 一 个 高 信 噪 比 的 脉冲 轮廓 , 将 其 作为 标准 轮廓 , 再 将 标准 轮廓 与 每 个 平均 脉冲 轮廓 做 互相 
关 , 得 到 脉冲 到 达 台 站 的 时 间 。 为 消除 地 球 运动 的 影响 ,通常 需要 将 太阳 系 质心 看 做 惯性 系 ， 
把 PSRCHIVE 得 到 的 所 有 到 达 时 间 转 换 为 到 达 太 阳 系 质心 时 间 ， 此 过 程 可 使 用 TEMPO2 实 
现 。 使 用 TEMPO2 将 本 次 数据 中 实际 观测 到 达 时 间 与 脉冲 星 计 时 模型 预 佑 到 达 时 间 进 行 拟 
合 ， 得 到 计时 残 差 。 转 换 时 用 到 的 太阳 系 星 历 为 DE405， 时 间 系 统 为 质心 坐标 时 TDB。 脉 
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冲 星 的 脉冲 相位 及 指数 恢复 过 程 分 别 如 式 CIO 和 式 (2) 所 示 ; 
pCt) = Go + vt +r Dt? + 2 jt? (D 
v(t) = vo(t) + Av, + Avst + = Ait? + Avge t/a oS 


X COD "B, pož r-0 MAINA: v. v. Bod SISSY BOR ATS BR = Bf 
导数 。 式 (2) P, Avy, Ady, A’, AIA KARA, Aa E Br SR 
未 衰减 的 值 ，Ava 为 恢复 幅度 ，zra 为 恢复 时 标 。 
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图 1 脉冲 星 计 时 数据 处 理 流程 


Fig.l Processing flow of pulsar timing data 


3 数据 处 理 结果 及 分 析 
3.1 预 处 理 结果 


HRT 观测 Crab 脉冲 星 的 某 次 数据 经 PSRCHIVE 处 理 ， 所 得 轮廓 图 如 图 
信 噪 比 约 为 33.445， 积 分 时 间 约 为 95 分 钟 。 在 此 轮廓 图 中 ， 可 清楚 看 到 Crab 脉冲 星 轮 廓 的 
三 个 主 成 分 : 分 别 位 于 相位 0.3 和 0.7 附近 的 主 脉 冲 和 中 介 脉 冲 ， 以 及 位 于 主 脉 冲 前 约 0.1 
处 的 前 兆 脉冲 。PSRCHIVE 处 理 所 得 Crab 脉冲 星 发 生 周 期 跃 变 前 的 自转 参数 如 表 1 所 示 。 
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图 2 Crab 脉冲 星 轮廓 图 
Fig.2 Profile of Crab pulsar 


表 1 周期 跃 变 前 的 自转 参数 
Tab.1 The parameters of pre-glitch 


Parameter Value 


| | 
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MJD 
R.A. (J2000) 
Dec. (J2000) 

v (s) 
ù (s?) 

ö (s?) 
DM (cm”°pc) 
UNITS 
CLK 
EPHEM 


3.2 计时 观测 结果 及 分 析 


频率 的 二 阶 导 数 妆 。 结 果 表 明 该 星 在 2019 年 7 
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58676 
05:34:31.921 (5) 
+22:00:52.16 (6) 
29.6169568869 (5) 
-3.683443 (8) x10? 
12611 (D x10” 

56.7983 (1) 
TDB 
TI(TAD 
DE405 


使 用 TEMPO2 拟 合 研究 数据 中 Crab FERA RME, B POSUERIT 7 SAB E 


月 23 日 (MDJ58687) 附近 发 生 了 一 次 周期 


跃 变 事件 ,计时 残 差 如 图 3 所 示 。 脉 冲 到 达 时 间 拟 合 得 到 自转 频率 及 其 一 阶 导 数 随时 间 的 变 
化 ,如 图 4 所 示 ， 其 中 上 面板 为 跃 变 后 的 自转 频率 与 其 平均 值 的 差 值 ， 下 面板 为 自转 频率 的 


一 阶 导数 随时 间 的 变化 。 拟 合 所 得 该 星 跃 变 后 的 参数 如 表 2 所 示 。 


本 次 跃 变 幅度 为 Avs/v = 17.9 (1) x10°, 


是 该 星 较 大 的 一 次 跃 变 事件 。 其 自转 增 量 为 


Av, = 5.33 (4) x10" Hz; 频率 一 阶 导数 的 变化 幅度 为 Avo =3.43 (4) x10; 一 阶 导数 变 
化 量 为 A 了 = -1.26 COD x10”s”。 对 周期 跃 变 后 的 分 析 表 明 此 次 周期 跃 变 伴随 着 一 次 指数 
恢复 过 程 ， 恢 复 幅 度 为 Aya= 4.69256078 (3) x10 7”， 恢 复 系数 O=Avg/Av, ~ 0.88， 恢 复 时 间 


Tq 283 (3) d. 


将 数据 处 理 结果 与 脉冲 星 周期 跃 变 统计 表征, 结果 进行 对 比 , 统计 表 中 记录 本 次 
跃 变 事件 大 小 为 Av, [v = 36.0 CD x10?， 两 者 跃 变 大 小 在 一 个 量 级 ， 数 值 上 的 差别 或 为 


HRT 观测 数据 在 跃 变 后 的 部 分 缺失 引起 。 
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图 3 ”计时 残 差 


Fig.3 Timing residuals 
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图 4 自转 频率 及 其 导数 随时 间 的 演化 


Fig.4 The evolution of rotation frequency and its derivative with time 


表 2 跃 变 后 的 参数 


Tab.2 The parameters of post-glitch 


Parameter Value 
MJD 58742 
R.A. (J2000) 05:34:31.921 (5) 


Dec. (J2000) 


+22:00:52.16 (6) 


v GO) 29.616957459 (1) 
v (s?) -3.68474 (6) x10? 
B (s?) 22 (D x10? 
DM (em?pc) 56.7983 (1) 
UNITS TDB 
CLK TT (TAD 
EPHEM DE405 


3 结 F 


中 国 科学 院 国 家 授时 中 心 吴 平 观测 站 40 米 射电 望远镜 对 Crab 脉冲 星 进行 长 期 计时 观测 ， 
监测 到 该 星 于 2019 年 7 月 23 日 附近 发 生 了 一 次 周期 跃 变现 象 。 通 过 处 理 及 拟 合 观测 数据 ， 
分 析 了 此 次 周期 跃 变 的 过 程 ， 结 果 表 明 该 星 跃 变 幅度 为 Avwv = 17.9 (1) x10”， 并 伴随 着 7 
4=8.3 (3) d 的 指数 恢复 过 程 ， 恢 复 系 数 O ~ 0.88。 目 前 ， 周 期 跃 变 的 物理 机 制 未 被 合理 解 
释 ， 脉 冲 星 计 时 观测 对 研究 脉冲 星 内 部 结构 具有 重要 意义 。 
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Abstract: Pulsar glitch, is a rare phenomenon, and is a probe for studying the internal 
structure of pulsars. From February to December 2019, time-of-arrival measurements have been 
collected for Crab pulsar with the 40-meter radio telescope at Haoping Observatory of the 
National Timing Center during pulsar timing observations. These measurements are analyzed by 
“timing solutions” with TEMPO2 fitting program. We have observed a glitch in the Crab pulsar. 
The glitch occurred around July 23, 2019 (MJD58687) when Crab pulsar underwent an increase in 
the rotation rate of Av,= 5.33(4)x107 Hz, corresponding to fractional size of Av,/v = 
17.9(1)x10?. This event was followed by an exponential recovery process with a degree of 
recovery Q ~ 0.88. The monitoring program of the Crab pulsar glitch, indicate that the 40-meter 
radio telescope of the National Time Service Center has a pretty good receiving sensitivity and 
system stability for pulsar timing observations, and contribute samples to study the mechanism of 
pulsar glitch. 
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